arXiv:1505.06487vl [astro-ph.SR] 24 May 2015 


A new sdO+dM binary with extreme eclipses and reflection effect 


A. Derekas 1,2 , P. Nemeth 3 , J. Southworth 4 5 6 , T. Borkovits 5,1 , K. Sarneczky 2 , A. Pal 2 , B. 
Csak 1 , D. Garcia-Alvarez 6,7,8 , P. F. L. Maxted 4 , L. L. Kiss 2,9 , K. Vida 2 , Gy. M. Szabo 1,2 , 

L. Kriskovics 2 


ABSTRACT 

We report the discovery of a new totally-eclipsing binary (RA=06 h 40 m 29 s ll; 
Dec=+38°56'52"2; J=2000.0; R max =17.2 mag) with an sdO primary and a strongly irra¬ 
diated red dwarf companion. It has an orbital period of P O rb=0. 187284394(11) d and an 
optical eclipse depth in excess of 5 magnitudes. We obtained two low-resolution classification 
spectra with GTC/OSIRIS and ten medium-resolution spectra with WHT/ISIS to constrain 
the properties of the binary members. The spectra are dominated by H Balmer and Hell 
absorption lines from the sdO star, and phase-dependent emission lines from the irradiated 
companion. A combined spectroscopic and light curve analysis implies a hot subdwarf tem¬ 
perature of T e ff(spec) = 55 000 ± 3000K, surface gravity of logg(phot) = 6.2 ± 0.04 (cgs) and 
a He abundance of log(nHe/nH) = —2.24 ± 0.40. The hot sdO star irradiates the red-dwarf 
companion, heating its substcllar point to about 22 500 K. Surface parameters for the companion 
are difficult to constrain from the currently available data: the most remarkable features are the 
strong H Balmer and C Ii-m lines in emission. Radial velocity estimates are consistent with the 
sdO+dM classification. The photometric data do not show any indication of sdO pulsations with 
amplitudes greater than 7 mmag, and Ha-filter images do not provide evidence of the presence 
of a planetary nebula associated with the sdO star. 
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1. Introduction 


Hot subdwarf stars are located between the up¬ 
per main sequence and the white dwarf (WD) se¬ 
quence in the Hertzsprung-Russell diagram. They 
are evolved, core-hclium burning, low-mass stars 
(M ~ 0.5 Mpt ) with very thin hydrogen envelopes 
( Heberl 20091) . Among hot subdwarfs, sdO stars 
(T e ff > 38000 K ) represent a significantly smaller 
fraction than sdBs (T e s < 35000 K). Spectroscop¬ 
ically, sdO stars show a large variety: the two main 
group s are the H-rich (sdO ) and He-rich (He-sdO) 
stars. IStroeer et al.l ( 20071) showed that a strong 
correlation exists among surface temperature, He 
and C, N abundances in He-sdO stars. 

Canonical stellar evolution theory predicts that 
sdO stars evolve from sdB stars. Whilst bina- 
rity is quite frequent amon g sdBs, with a bi¬ 


nary fraction of about 50% ( Maxted et ahlliToOll 
Napiwotzki et al.1 2004 ). the fraction of binary 
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He-sdOs is very low ( Kawka et al.l I20T5I) . Many 
sdO binaries are associated with a planetary neb¬ 
ula (PN), like UU Sge , V477 Lyr and BE UMa 
dPollacco fe Belli 1994 Pollacco. Bell fe Hilditch 


1994 ; Afsar fc ibanoglu 2008h . The sdO stars in 
these binaries are hotter and more massive stars 
in the immediate post giant branch stage. Bi¬ 
naries with compact sdO stars evolve from sdB 
binaries and the ~120 Myr sdB lifetime is long 
enough for their PNe to attenuate and become 
hardly detectable ( Aller et al. 20151) . 

lHan et ah ( 2002 . 12003 ) performed binary pop¬ 
ulation syntheses and identified several evolution¬ 
ary channels that lead to the formation of hot 
subdwarf stars. In close binaries that evolve 
through one or more common envelope (CE) 
phases ( Paczvnskill976l) . the secondary is engulfed 
by the atmosphere of the primary while it is on 
the red giant branch. 

As the secondary spirals inward due to tidal 
friction, the red giant loses mass. By the end of 
the CE phase the primary loses most of its enve¬ 
lope and the binary orbital period shrinks to a few 
hours. If the core gained enough mass for He ig¬ 
nition during the preceding evolution the primary 
experiences a He-flash and settles on the extreme 
horizontal branch. In case the core mass is in¬ 
sufficient for He burning, the primary evolves as 
a low-mass pre-white dwarf. In both cases the 
common envelope is ejected during the final stage 
of CE evolution and a very close binary remains 
( Taam fe Ricker 20061) . While this theory can ex¬ 
plain the mass loss required for the formation of 
hot subdwarfs, it also needs a precise timing be¬ 
tween mass loss and the core helium flash. Eclips¬ 
ing hot subdwarf binaries with irradiated compan¬ 
ions can give insight into the details of these pro¬ 
cesses, making such binaries fundamental to un¬ 
derstand CE evolution. 


Here we report the discovery of a new eclips¬ 
ing binary with an sdO primary and a strongly 
irradiated red dwarf companion. After describ¬ 
ing the observations, we discuss the spectral mod¬ 
elling and determination of the stellar and orbital 
parameters from the spectroscopic and light curve 
analysis. 


2. Observations 


Konkoly J064029.1+385652.2 (RA=06 h 40 m 29 s ll 
Dec=+38°56 , 52"2; J=2000.0; R max =17.2 mag; 
hereafter J0640+3856) was discovered serendipi- 
tously during regular astrometric observations of 
minor planets. On one of the images the object 
completely disappeared, suggesting a sudden deep 
eclipse. We started to monitor J0640+3856 using 
the 0.6/0.9/1.8 m Schmidt telescope at Piszkesteto 
Observatory. We took CCD photometric obser¬ 
vations on seven nights between December 2013 
and February 2014 using Johnson/Bessell V and 
Cousins R c filters, and also without filters. The 
telescope was equipped with an Apogee ALTA-U 
4kx4k CCD camera. The observations revealed 
that the period is 0.187 d and there is a strong re¬ 
flection effect with an amplitude of ^0.5 mag. The 
sudden and deep primary minimum suggested that 
the system might contain a hot and small primary 
star eclipsed by a cool secondary object. 


We obtained further observations via the ser¬ 
vice program on the 4.2 m William Herschel Tele¬ 
scope (WHT) on La Palma. In February and 
March 2014 we observed two primary eclipses 
with Sloan r' and %' filters and a secondary 


eclipse in i' band, using the ACAM imager 
( Benn. Dee fc Agdca 2008 0. The data were re¬ 
duced with standard procedures. The flat minima 
indicated that the eclipse depths are 6 and 5 mag, 
respectively, and the totality of the primary eclipse 
lasts for about 5.7 min, while the eclipse duration 
is approximately 24 min. 


On March 1, 2014, we obtained two spectra us¬ 
ing the 10.4m Gran Telescopio Canarias (GTC) 
and OSIRIS spectrograplQ on La Palma. Due to 
the highly variable seeing we used a 1.23 arcsec 
slit providing a dispersion of AA « 1.6 A/pixel 
and AA ~ 4.2 A/pixel with the R2500V and the 
R1000B grisms, respectively. These low- and 
medium-resolution spectra were taken between 
3850-7400 A and 4400-6000 A, and reached a 
signal-to-noise ratio (SNR) ~ 30 with 540 s expo¬ 
sure times. 


We acquired phase-resolved spectroscopy on 
March 6, 2014, using the WHT and the dual-beam 
ISIS spectrograpto This was operated with the 


1 http: //wrn. gtc.iac.es/instruments / osiris 

2 http: //www.ing.iac.es/Astronomy/instruments/isis 
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R600B/R600R gratings and a 1 arcsec slit, pro¬ 
viding a resolution of A A = 1 A/pixel in the blue 
and covering the 3800-5200 A and 6200-7800 A re¬ 
gions. Since the Ha line was weak, we did not use 
the red region in the analysis. We took ten spec¬ 
tra between orbital phases </? = 0.47 and 0.82 with 
600 s exposure times. The average SNR of these 
spectra is 25. 

All spectroscopic data were reduced with our 
IRAF based data reduction pipeline. Bias and 
flat field corrections were done with the CCD- 
PROC task, raw spectra and arc calibrations were 
extracted with APALL and wavelength calibrated 
with the identify/reidentify tasks. We iden¬ 
tified —150 CuArNe lines in the WHT calibra¬ 
tion data, that allowed for an eighth order Leg¬ 
endre polynomial dispersion function with a root- 
mean-square (RMS) —0.1 A. The GTC HgArNeXe 
calibrations allowed us to identify —50-100 lines 
and use a sixth order dispersion function with 
RMS-0.6 A. 

In April 2014 we obtained 1.5 h CCD photome¬ 
try with the 1 m RCC telescope at the Piszkesteto 
Observatory, using an FLI camera with a field of 
view of 9.4'x9.4', with a Sloan r' filter and 2x2 
binning giving a plate scale of 0.27"/pixel. In ad¬ 
dition, we obtained — 21 h of fast photometry (5 s 
exposure times) in December 2014 and February 
2015 using the same telescope with the OCELOT 
camera (an Andor iXon+888 EMCDD camera) 
without filter, in order to search for any oscilla¬ 
tions from the sdO star. We did not detect any 
significant periodic signal with amplitudes greater 
than 7 mmag. 

In order to detect whether any PN could be as¬ 
sociated with the sdO, we took 27 Ha-filter images 
with 120 s exposure times with the same telescope 
in April 2014. The composite image did not show 
any sign of a PN. 


3. The orbital period 


We checked the stability of the orbital pe¬ 
riod. In order to determine accurate times of 
minima the WHT light curves of the two pri¬ 
mary and one secondary minima were used as tem¬ 
plates, and these template curves were fitted to 
th e other, less accurate light curves as described 


Borkovits et al. ( 2015h in details. Thus, we ob¬ 


tained 13 primary and 14 secondary times of min¬ 


ima, which are listed in Table Q] and plotted in 
Fig. [T] For the period variation analysis, however, 
only the primary minima were used, with the one 
exception of the outlying point at JD 2 456691. 
Therefore, we used twelve primary minima ob¬ 
tained between JD 2 456 684 and JD 2 457 068 to 
calculate the following ephcmeris: 

T min , i = 2 456 709.58464(2)+0 < ?187284762(81)xE, 

where the epoch is given in BJD (TDB). We did 
not find any evidence for a variation in the period, 
however, it should be noted that the scattering 
of the individual points exceeds significantly the 
estimated statistical errors. 


4. Spectroscopy 


4.1. Modelling the composite spectrum 

The first light curve of J0640+3856 clearly sug¬ 
gested a HW Vir type eclipsing binary with a 
strongly irradiated cool companion. In such bi¬ 
naries the orbit is small and the companion is rel¬ 
atively large to produce a prominent reflection ef¬ 
fect in the light curve and provide full eclipses of 
the primary. 

Our classification spectra showed strong emis¬ 
sion lines of C superimposed on the shallow ab¬ 
sorption spectrum of a hot star. J0640+3856 
is therefore a double lined spectroscopic binary, 
making it a target with high potential. Disentan¬ 
gling the spectrum of the hot spot usually requires 
high-contrast high-SNR spectroscopy and a chal¬ 
lenging analysis procedure, similar to direct ex¬ 
oplanet spectroscopy. The spectra of irradiated 
companions have been disentangled and analysed 
in several hot s ubdwarf binaries, e.g. EC 11575 
and V664 Cas (lExter etldl l2005lb WD0137-349 


( Maxted et al. 20061) andAA Dor (Vuckovic et al 
20081) . However, these analyses were limited to 
relative line strength and radial velocity measure¬ 
ments based on Gaussian line profiles. Although 
grids of irradiated M-dw arf spectra are available 
(e.g. iBarman et al. 2004 ). their applicatio n to fit 
observed data is limited. Wawrzvn et al.1 ( 20091) 
presented a self-consistent model atmosphere anal¬ 
ysis for the hot subdwarf binary UU Sge. Here we 
foll ow their met hods based on the description in 
Gunther fe Wawrzvnl ( 201ll) . 

J0640+3856 showed a double lined composite 
spectrum immediately in low-resolution moderate- 
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Table 1: Times of minima of J0640+3856. 


Primary minima 

Secondary minima 

BJD (TDB) 

-2 400 000 

Cycle 

no. 

std. dev. 
(d) 

instrument 

BJD (TDB) 

-2 400 000 

Cycle 

no. 

std. dev. 
(d) 

instrument 

56684.301289 

-135.0 

0.000004 

Schmidt Rc 

56654.616610 

-293.5 

0.000019 

Schmidt unfilt. 

56684.488447 

-134.0 

0.000004 

Schmidt Rc 

56657.424647 

-278.5 

0.000019 

Schmidt unfilt. 

56691.417455* 

-97.0 

0.000031 

Schmidt V 

56684.394740 

-134.5 

0.000019 

Schmidt Rc 

56692.354351 

-92.0 

0.000004 

Schmidt V 

56689.265080 

-108.5 

0.000032 

Schmidt V 

56693.290882 

-87.0 

0.000004 

Schmidt Rc 

56691.325691 

-97.5 

0.000007 

Schmidt V 

56695.350946 

-76.0 

0.000001 

WHT r’ 

56691.509921 

-96.5 

0.000017 

Schmidt V 

56709.584600 

0.0 

0.000002 

WHT i’ 

56692.260600 

-92.5 

0.000019 

Schmidt V 

56775.321665 

351.0 

0.000021 

RCC r’ 

56692.448527 

-91.5 

0.000024 

Schmidt V 

57022.537442 

1671.0 

0.000005 

RCC unfilt. 

56693.385087 

-86.5 

0.000023 

Schmidt Rc 

57067.298585 

1910.0 

0.000002 

RCC unfilt. 

56729.530347 

106.5 

0.000016 

WHT i’ 

57067.485826 

1911.0 

0.000003 

RCC unfilt. 

57022.631307 

1671.5 

0.000021 

RCC unfilt. 

57068.234937 

1915.0 

0.000003 

RCC unfilt. 

57067.391972 

1910.5 

0.000015 

RCC unfilt. 

57068.422263 

1916.0 

0.000002 

RCC unfilt. 

57068.328437 

1915.5 

0.000013 

RCC unfllt. 





57068.515991 

1916.5 

0.000018 

RCC unfllt. 


Note. Secondary minima, and the one primary minima denoted by * were not included in the period study. 
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Fig. 1.— Eclipse Timing Variation diagram of J0640+3856. Red circles represent the primary and blue 
boxes the secondary minima. For a better visibility, secondary minima with a large scatter are not shown 
here. The black line represents the linear fit for the primary minima, which we used to calculate the ephemeris 
and orbital period given both in the text and the figure. It shows that the period is well determined and 
constant over one year. 
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T 0 = 2456709.58464(2) BJD (TDB) P = 0.187284762(81 ) c 
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Fig. 2.— Spectral decomposition of the WHT/ISIS spectrum at ip = 0.647. The spectral features can 
be described with the superposition of a 55 000 K sdO star and an irradiated dM companion heated to 
~22 500 K. 


SNR spectroscopy, as achievable with 4m-class 
telescopes, despite being fainter by ^6 magnitudes 
in V than AA Dor. This suggests that both com¬ 
ponents are more extreme: the sdO star must be 
hotter and the companion must be a larger, prob¬ 
ably earlier-M star that suffers a stronger irra¬ 
diation. The most remarkable spectral features 
are the distorted, phase-dependent H Balmer line 
profiles, the strong He n absorption indicating an 
sdO star and C ii-m lines in emission from the ir¬ 
radiated companion. These features defined the 
starting models in our analysis. The WHT ob¬ 
servations include the secondary minimum, what 
makes them ideal to investigate the irradiated 
hemisphere of the companion. 


Our spectral analysis was based on the steepest- 
descent itera tive binary fit proce dure, XTgrid, 
developed by Nemeth et al. (I2012h . XTgrid em¬ 


1995; 


ploys Tlusty/Synspe c (jHubenv fe Lanz 
Lanz fc Hubcnvl l2003h non-LTE model atmo¬ 


spheres and synthetic spectra to reproduce com¬ 
posite binary spectra. The strong irradiation al¬ 
lowed us to model both components in J0640+3856 
with Tlusty. The model atmospheres include 
H, He, C, N and O opacities consistently in the 
atmosphere structure and synthetic spectrum cal¬ 
culations for both stars. Although we could not 
identify CNO absorption lines in the sdO star and 
kept their abundances fixed at log(nX/nH) = —6, 
we included CNO opacities, because they have an 
effect on the tempe rature structu re of non-LTE 
atmosphere models (IWerncrlll996lf . 


We started an XTgrid process for the sdO star 
and another for an M dwarf that is irradiated by 
the sdO. These processes iteratively updated the 
stellar parameters and fit the observed composite 
spectrum together. As a first approximation we 
applied an isotropically irradiated model for the 
inner hemisphere of the companion. We assumed 
a black-body energy distribution of the irradiat¬ 
ing flux with appropriate temperature for the sdO . 
The geometric dilution factor (W) is (IMihalas 


197811 : 




where R 2 is the radius of the secondary and a is 
the semi-major axis. The dilution factor describes 
the strength of irradiation. We used an empiri¬ 
cal dilution factor W e = 0.05 that was found for 
the substellar point in the similar system AA Dor. 
The empirical dilution factor is a fudge factor in 
our analysis and substantially larger than the ge¬ 
ometrical dilution factor. The difference is prob¬ 
ably due to the non-planckian sdO spectrum and 
due to the fact that the spectrum of the irradi¬ 
ated companion cannot be fully described with the 
conditions i n the subs tellar point, as suggested by 
Gunther & Wawrzvn ( 201lh . 


The combined spectrum at ip = 0.647 has a 
SNR of ~ 45, therefore we modelled this spec¬ 
trum as shown in Fig. [2] The flux contribution 
of the hot spot was determined for this spectrum 
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Fig. 3.— Top: GTC/OSIRIS spectra of J0640+3856 (grey) with best fit models (black). Bottom: Best fit 
models (black) for the phase-dependent radial velocity corrected composite spectra obtained with WHT/ISIS 
(grey). We applied a spectral decomposition of the middle spectrum at orbital phase = 0.647. All the 
other models were calculated from this decomposition by shifting the components in radial velocity and 
scaling them by the orbital phase to account for the visibility of the substellar point of the companion. 
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and we approximated it for the other spectra by 
scaling with orbital phase according to Fdu/Fsd ~ 
0.43 sin 2 (7 rip). The Tlusty model provides the 
temperature structure of the irradiated compan¬ 
ion. The temperature decreases steeply inward, 
reaches a minimum and progressively increases 
again. The minimum temperature is T = 22 500 K 
at optical depth r « 0.6, which we associate with 
the photospheric effective temperature of the sub- 
stellar point. The parameters derived from the 
spectroscopic modelling are listed in Tabic [2] and 
suggest an sdO+dM binary. 

Next we shifted the components in velocity 
space to reproduce the phase-dependent compos¬ 
ite spectra. Fig. [3] shows the best fits for the ten 
WHT blue arm observations. These fits show that 
all H line profiles are contaminated by the emis¬ 
sion lines from the companion and the relative 
strength of this contamination increases towards 
the Balmer series limit. The inverted temperature 
structure of the companion results in high excita¬ 
tion lines compared to the photospheric tempera¬ 
ture. Such emission lines of C ill and Si ill have 
been observed from the irradiated companion in 
other systems (e.g.: V477 Lyr, AA Dor). There¬ 
fore, a He 11 4686 A emission line may be expected 
as well. However, this line forms between very 
high lying levels in the He ion and at the low He 
abundance of J0640+3856 we did not find any con¬ 
tribution. The fact that the He II 4686 A line forms 
exclusively in the sdO star allowed us to measure 
its radial velocity. The CII 4267 A line, the C ill 
blend near 4650 A and 4070 A are in emission and 
come from the irradiated companion. Based on a 
selection of these lines we could also estimate the 
radial velocity of the companion. 

4.2. Radial velocity and stellar masses 

Even though precision radial velocity measure¬ 
ments would require high-dispersion and higher- 
SNR spectra, we attempted to estimate radial ve¬ 
locities from the WHT/ISIS spectra, because the 
amplitudes put a useful constraint on the mass 
ratio. First, we checked the dispersion correction, 
which was based on 150 lines of the CuAr and 
CuNe calibration lamps, and tested the stability 
on the CailK line. This line forms mostly in the 
interstellar medium and therefore its radial veloc¬ 
ity can be used to check dispersion offsets. We 
found the K line to be consistent in all our spec¬ 


tra, but systematically offset by 40kms _1 . Most 
of this offset can be explained with the barycen- 
tric velocity correction vbc = —26.5 kms -1 dur¬ 
ing the observations. The interstellar extinction 


towards J0640+3856 is E (B - V) = 0.124 ± 0.005 
( Schlaflv fe Finkbeiner 201lh . therefore we asso¬ 
ciate the remaining radial velocity shift of the 
CailK line with the relative velocity of the in¬ 
terstellar material with respect to the Solar sys¬ 
tem. We performed two independent experiments 
to measure velocities and limited our measure¬ 
ments to the He 11 4686 A line for the sdO and to 
the C ill blend between 4647-4652 A for the dM 
companion as these are the strongest undistorted 
features in the spectra. 


In the first method we combined each consecu¬ 
tive spectrum, like those at ip = 0.622 and 0.673 
to get the one at 0.647 in Fig . OH and applie d 
a Savitzky-Golay filter ( Savitzkv fc Golav 1964 ). 
Then we measured the SNR of each of these com¬ 
bined spectra. Next, we determined the radial ve¬ 
locity with an iterative chi-square minimization 
cross-correlation procedure. We decided to use 
the preceding spectrum as a self-template to re¬ 
duce systematic effects in the chi-square method. 
For each observed spectrum we performed 100 ra¬ 
dial velocity measurements to obtain a mean value 
and standard deviation. We resampled the spec¬ 
tra according to the SNR before each individual 
measurement. 


The second method was based on visual inspec¬ 
tion. This allowed us to compare the observations 
to synthetic spectra directly and derive absolute 
radial velocities. 


We found a good match between the two meth¬ 
ods for the sdO star. However, for the dM com¬ 
panion we found that the cross-correlation method 
underestimated the radial velocity by about 40%. 
The most plausible reason for this is orbital smear¬ 
ing that makes the spectral lines considerably shal¬ 
lower and broader, while the decreasing contribu¬ 
tion of the companion also makes the last data 
points less reliable. Similarly for the sdO star, 
the last data points show a decreasing radial ve¬ 
locity in conflict with the photometric orbital pe¬ 
riod. We attributed these inconsistencies to a 
lower data quality and disregarded the last three 
spectra in the radial velocity measurement. The 
symmetric appearance of the reflection effect and 
secondary minimum suggest a circular orbit and 
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a co-rotating companion. Assuming synchronised 
rotation we found a projected rotational velocity 
^rotsini = 49.2 km s -1 using R 2 = 0.2 R 0 from 
the light curve solution (see Sect. ED). As the centre 
of mass of the companion does not coincide with 
its centre of light we needed to correct the radial 
velocity with the projected rotational velocity at 
the centre of light. We calculated this correction 
as — 29kms -1 at <p = 0.75. 



Fig. 4.— Radial velocity curves of the components 
of J0640+3856. 

The six data points in Fig. |4] show a clear 
trend and suggest the corrected radial velocity 
semi-amplitudes: K\ = 80 ± lOkrns -1 , K 2 = 
257 ±25kms _1 and a barycentric system velocity 
Vj = 61 ± 8 km s -1 . From the radial velocity semi¬ 
amplitudes we found a mass ratio of q = M 2 /M 1 = 
Ki/K 2 = 0.31 ± 0.04. Assuming a circular orbit, 
the radial velocity semi-amplitudes define the pro¬ 
jected semi-major axis asini rs ( K\ + K 2 ) P/2tt, 
where P is the orbital period and i is the incli¬ 
nation, which can be obtained from the eclipsing 
light curve solution. The observed orbital period 
and this semi-major axis define the total mass by 
Kepler’s third law. 

Using the inclination from the light curve solu¬ 
tion i = 87.11° and the radial velocity amplitudes 
the semi-major axis is a = 0.006 AU (869 000 km). 
Then the total mass from Kepler’s third law is 
M tot = 0.744 M 0 . 


According to our models the strong irradiation 
transforms the secondary atmosphere completely. 
The irradiated models are insensitive for the un¬ 
perturbed temperature of the companion, there¬ 
fore we cannot assess the night side temperature 
from the currently available irradiated spectrum. 
In turn, the models show that the Balmer emission 
lines are sensitive for the surface gravity (pressure) 
and to reproduce the line profile variations they 
require a log <7 >4.5. Assuming that the mass of 
the subdwarf is close to the canonical 0.5 M 0 , the 
measured surface gravity log g = 4.9 of the com¬ 
panion is most consistent with a mid-M type star. 


5. Light curve analysis 


We carried out a three-band simultaneous 
light curve analysis with the recent l y devel ¬ 
oped Light curvefactory code ( Borkovits et al.l 
I2013LI 2 OI 4 I I. We chose the Sloan i' and r' light 
curves obtained using WHT/ACAM, which con¬ 
tain two flat primary minima offering strong ge¬ 
ometrical constraints not only upon the incli¬ 
nation and the relative radii of the stars, but 
also, via Kepler’s third law, the surface gravity of 
the stars given an approximate estimate of their 
mass. The third light curve used was the Rc~ 
band dataset obtained with the Schmidt telescope 
at the Piszkesteto Observatory, which covers more 
than a full orbit of the system, including two pri¬ 
mary minima. Despite a lower quality and the 
absence of measurements during the deepest parts 
of the primary minima, these data significantly 
improved the fit of the reflection effect and also 
contain useful information on the ellipsoidal vari¬ 
ation or its absence that can be used to further 
constrain the mass ratio. 


There is no indication for eccentricity from 
the spectroscopic orbit or the phase of secondary 
eclipse so we assumed that the orbit is circular, 
as expected for such a short-period binary star. 
Assuming a spherical primary, and a marginally 
oblated secondary star (which latter assumption 
was found to be reliable from a preliminary analy¬ 
sis of the eclipse geometry), it is a suitable approx¬ 
imation that the fractional radii of the components 
relate directly to the observable quantities of the 
full and totality-phase durations of the primary 
occultations. In such a way the whole system ge¬ 
ometry was determined except for one free param- 












-0.20 0.00 0.20 0.40 0.60 

Phase 


Fig. 5.— Observed light curves and solutions with their residuals for the Sloan i', r’ (WHT) and Rc 
(Piszkesteto Observatory) passband measurements. 


eter, namely i. Therefore we were able to express 
the adjustable dynamical and geometrical param¬ 
eters using the binary period (P), epoch (To), the 
orbital inclination (/) and the mass ratio ( q ). Al¬ 
though the mass ratio is also known from the ra¬ 
dial velocity solution (see Sect. 14.21) . we made fits 
both with it freely adjusted and fixed to the spec¬ 
troscopic value, as a consistency check. Due to the 
detached configuration and nearly spherical stellar 
shapes, the mass ratio has only a minor influence 
on our light curve solution. 

Considering the atmospheric properties, the ef¬ 
fective temperature of the sdO primary was fixed 
to the value obtained from the spectroscopic anal¬ 
ysis (T e ffi = 55 000K), while T e ff 2 ,d (night side of 
the secondary component) was the fifth adjusted 
parameter. The other atmospheric parameters - 
limb darkening, gravity brightening coefficients, 
bolometric albedos and abundances of the sdO 
primary - were also kept fixed. For limb darken¬ 
ing we applied the logarithmic law, and the coeffi¬ 
cients were calculated according to the passband- 


depen dent precomputed tables of t he phoebe 
team ( Prsa fc Zwittcrl[2005l Prsa[2011 ) which were 
based on the tables of Castelli fc Kurucd ( 20041 ). 

On the other hand, most of the atmospheric 
parameters of the secondary star were involved 
into the fitting process. This was done because of 
the high irradiation which results in such a high 
temperature as T = 22 500 K at the substellar 
point that there are significant d eviations from 


LTE models as discussed e.g. in iBarman et al 
( 20041) . Amongst other consequences, this may 
change the full atmosphere to a radiative one, 
as was o bserved for other si milar systems (e.g. 
AA Dor; iHilditch et al.ll2000l ). and may even re¬ 
sult in limb-brightening ( i.e. negative limb dark - 
ening coefficients, see e.g. IPollacco fe Bell ( 1994 ). 
for V477 Lyr). 


Finally, the passband luminosities of the pri¬ 
mary component were also re-calculated in each 
step. A x 2 minimization was carried out us¬ 
ing firstly the Levenberg-Marquardt differential 
corrections algorithm, and then a final refine- 


3 Downloaded from the site http://phoebe-project.Org/l.0/ 
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ment with a grid-search method. The resulting 
light curve solutions, and their residual curves are 
shown in Fig. [5] The values of the free param¬ 
eters found from the best fit to the light curves 
are given in Tabic [2J together with the fixed pa¬ 
rameters in the least-squares fit and some useful 
derived parameters. 

The photometric mass ratio was found to be 
significantly smaller (q = 0.19 ± 0.02) than the 
spectroscopic one. On the other hand, compar¬ 
ing the two solutions (with adjusted or fixed mass 
ratios), we found that all the other adjusted pa¬ 
rameters, and also the y 2 have remained within 
their lcr values. This is a consequence of the small 
amplitude ellipsoidal variation, due to the minor 
oblateness of the secondary star and, therefore, 
the weak dependence of the light curve on q. We 
therefore decided to keep the solution obtained 
with the fixed spectroscopic mass ratio and the de¬ 
rived astrophysical quantities were computed ac¬ 
cordingly. Considering the other quantities com¬ 
mon to the spectroscopic and photometric analy¬ 
ses, while the light curve solution clearly confirms 
the spectroscopic temperature of the secondary, 
the local gravities were found to be higher. 

Using the spectroscopic mass ratio we found 
a mass of 0.57 ± 0.14 M 0 for the sdO star and 
0.18 ± 0.05 Mq for the companion, in agreement 
with a mid-M dwarf (M6V) classification. We 
note that although the spectroscopic mass ratio 
is higher, yet consistent with our light curve anal¬ 
ysis, the poor quality of our spectra means it is 
not precisely determined. 

Comparing the spectroscopic and photometric 
logg values of the sdO primary, despite that the 
lower spectroscopic value is closer to expectations, 
the photometric result is evidently the more ro¬ 
bust. This is because it is determined purely by 
the system geometry via the eclipse durations, 
which give the relative radii as a the function of 
i. Then, by the use of Kepler’s third law, one can 
see that g pr ; ~ + g) 2 / 3 . Therefore, even 

a 100% error in the mass of the primary would 
result in a 0.1 dex discrepancy in its log g value. 
This is valid only for spherical stars with negligible 
tidal and rotational effects, but our solution (i.e. 
the low values of the fractional radii) is consistent 
with the assumption that these effects play only a 
minor role in the system. As a consequence, we ac¬ 
cept the high log<? values obtained from the light 


curve analysis instead of their spectroscopic value 
and conclude that the sdO primary is a compact 
object on its way to the WD cooling sequence. 


Turning to the adjusted atmospheric parame¬ 
ters of the secondary, its bolometric albedo (A 2 ) 
was found to be greater than unity. This is not an 
unphysical solution, but implies that light from 
outside of a given photometric passband is be¬ 
ing reprocessed and re-emitted in the passband. 
A similar situation has bee n found before for 
other systems (e.g., ISouthworth et alj 2011 . for 
KIC 10661783, a totally eclipsing binary wi th a 
5 Scuti component and lAlmeida et al.112012L for 
the eclipsing sdOB+dM binary V1828 Aql). We 
also found an unusually high gravity darkening co¬ 
efficient (@ 2 ) for the secondary. The large uncer¬ 
tainty, however, makes this result ambigous. We 
note that the relatively poor fit to the data in 
secondary minimum is likely due to the simplis¬ 
tic treatment of the most highly-irradiated part 
of the atmosphere of the secondary component, 
which is eclipsed during secondary minimum. A 
complete physical description of the phsyics in this 
atmosphere is not available in our code, or in other 
commonly available light curve modelling codes. 


6. Summary and Conclusions 


We have discovered an sdO+M6V eclipsing 
binary (RA=06 h 40 m 29 s ll; Dec=+38 0 56'52"2; 
J=2000.0; R max =17.2 mag) that shows 6-mag 
deep primary eclipses and a ^0.5 mag reflection ef¬ 
fect. These are the most extreme variations among 
all HW Vir type binaries known to us and the pri¬ 
mary minimum is even deeper than that of NN Ser, 
a well- known white dwarf w ith an extreme eclipse 
depth (IHaefner et al . 20041 ). With photometric 
and spectroscopic follow-up we constrained the 
atmospheric properties of the components and the 
binary orbit. Although the specific spectral fea¬ 
tures and the effective temperature (T e m = 55 000 
K) classify the primary component as an sdO star, 
the surface gravity is at the upper limit of sdOs 
(log g = 6.2 cgs) and the radius (R = 0.096 Rq) is 
smaller than for normal sdO stars. These param¬ 
eters place the primary component in a special 
position, suggesting the sdO star is a pre-WD, 
similar to BE UMa which is classified as a border¬ 
line obj ect between sdO subd warfs and DAO white 
dwarfs ( Ferguson et al. 1999 ). The non-detection 
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Table 2: Stellar and orbital parameters derived from the spectroscopic and eclipsing light curve analysis. 

orbital parameters 
0.18728550 ±0.00000005 
2 456 709.584565 ± 0.000013 
1.24857 ± 0.09976 
0.0 

87.11 ±0.03 
0.31 ±0.05 

_ 0.19 ±0.02 _ 

stellar parameters 

Primary Secondary 

fractional radii d 

0.07646 ± 0.00053 0.15703 ± 0.00055 

0.07646 ± 0.00053 0.15838 ± 0.00055 

0.07646 ± 0.00053 0.16203 ± 0.00065 

0.07646 ± 0.00053 0.16130 ± 0.00063 


absolute stellar parameters 


M d (M 0 ) 

0.567 ±0.138 

0.177 ±0.051 

R d (R e ) 

0.0955 ±0.0077 

0.1985 ±0.0159 

T e ff (spec) (K) 

55 000 ± 3 000 

4000™ 

T eff (phot) (K) 

55 000 (±3 000) 

4648 ± 55 (±259) 

-^bolCL©) 

73.692 ± 11.819 (±19.955) 

0.016 ±0.004 (±0.007) 

logg (spec) (cgs) 

5.97 ±0.30 

4.9 ±0.5 

log g d (phot) (cgs) 

6.23 ±0.04 

5.11 ±0.07 


chemical abundances from spectroscopy 

log(nHe/nH) 

-2.24 ±0.4 

-2.4 > 

log(nC/nH) 

-6.0 

-2.0 > 

log(nN/nH) 

-6.0 

-5.0 > 

log(nO/nH) 

-6.0 

-3.5 > 


flux 

ratio from spectroscopy 

(FdM/Rsd ) max at 4500 A 


0.43 ±0.06 

atmospheric model-dependent parameters 

*^bol 

0.231 

0.300 ±0.01 

2 /bol 

0.148 

- 

A 

1.0 

1.09 ±0.02 

P 

1.0 

2.80 ±0.70 

Xi 

0.160 

0.47 ±0.02 

Vi 

0.108 

- 

Li/ ± ±2)(*) 

0.827 

0.173 

x r 

0.184 

0.50 

Vr 

0.122 

- 

Li/(Li ± ±2)(f) 

0.986 

0.014 

%Rc 

0.179 

0.50 ±0.02 

VRc 

0.119 

- 

L i /{L 1 + L 2 ){Rc) 

0.984 

0.016 


Notes. (1) Parameters without uncertainties were kept fixed, or adopted from precomputed tables. (2) Parameters 
subscripted with d are derived parameters. (3) Second uncertainties in parentheses were calculated setting the un¬ 
certainty of the (fixed) primary effective temperature to be 8T e m = 3 000K (i.e. its spectrocopic uncertainty). (4) 
x, y , A and /3 denote linear and logarithmic limb darkening coefficients, bolometric albedos and gravity brighten¬ 
ing exponents, respectively. (5) In the case of the passband-dependent fractional luminosities [Li/(L i ± ± 2 )], the 
reflection/irradiation effect was taken into account. 


fpole 
fside 
T" point 
Hiack 


Porb (days) 
Tmini (BJD) 
a d (Rq) 
e 

i(°) 

qd 

y spec 

<? P hot (unused) 
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of a PN around J0640+3856 also supports the 
evolved hot subdwarf (post-sdB) scenario. 

We have constructed a simple model to repro¬ 
duce the spectral contribution of the irradiated 
companion. Although this model is optimized for 
the substellar point, it represents the day side of 
the companion well, suggesting that the strong ir¬ 
radiation heats up the entire inner hemisphere ho¬ 
mogeneously. 


Our results suggests that the secondary com¬ 
ponent may be inflated by only a few per¬ 
cent; like i n the cases of similar close binaries 
( Afsar fc i bano glu 12008 1. We estimate that the 
substellar point of the red dwarf is heated to about 
22 500 K. The heat transport of these inflated stars 
is ineffective, so the large temperature difference 
between the day and night sides is preserved over 
long timescales ( Ritter. Zhang fc Kolbl 2000 (1. 

We conclude that the most probable compan¬ 
ion spectral type is mid-M. A later-type or more 
compact companion would be unable to reproduce 
the eclipses while an earlier type and more massive 
companion would be inconsistent with the radial 
velocity curve. 


The biggest advantage of J0640+3856 is that 
it is a double-lined spectroscopic binary. Spec¬ 
troscopic observations covering the full orbital cy¬ 
cle will yield more precise parameters (especially 
masses) for the components, as well as an opportu¬ 
nity to monitor and analyse the changing features 
in the spectra caused by the reflection effect. 


The J0640+3856 system is a good analogue to 
study interactions in planetary systems with hot 
Jupiters. Both the illumination effect in the pri¬ 
mary minimum and the thermal radiation and re¬ 
flected light disappearance and reappearance in 
the secondary minimum, is similar as is the lu¬ 
minosity ratio. 


The newly discovered J0640+3856 is a unique 
laboratory in several aspects and opens opportuni¬ 
ties to make further steps to understand the evolu¬ 
tionary history of post-common envelope binaries. 
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